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Abstract 

Physical processes ranging from the Lamb shift to the energy loss dE/dx of 
a charged particle traversing a plasma entail processes that occur over a wide 
range of energy or length scales. Different physical mechanisms dominate at 
one or the other end of this range. For example, in the energy loss problem, 
soft collisions that are screened by collective effects are important at large 
distances, while at short distances hard collisions are important where the 
exact details of the single-particle interactions must be taken into account. We 
introduce a novel application of dimensional continuation. The soft processes 
dominate at all scales when the spatial dimension u is less than 3, and we use 
them to compute the result to leading order for v <2>. On the other hand, the 
hard processes dominate at all scales for > 3, and we use them to compute 
the result to leading order for these spatial dimensions. We then explain why 
the sum of the analytic continuation of these disparate mechanisms yields the 
correct leading-order result for the physical limit at = 3 dimensions. After 
applying this new method to the energy loss problem in some detail, we then 
show how it also provides a very short and easy way to compute the Lamb 
shift. 
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I. INTRODUCTION 



The purpose of this paper is to introduce a new apphcation of dimensional continuation 
to physical problems that involve simultaneously both short-distance, ultraviolet processes 
and long-distance, infrared processes. Such problems typically involve a logarithm of a large 
ratio of two different scales. Although the overall coefficient of the logarithm is usually easy 
to compute, the constant under the logarithm is often difficult to find. Our new method 
makes the computation of this constant under the logarithm easy. We shall illustrate the 
simplicity of the method and its wide range of applicability by using it to solve two very 
diverse physical problems, the dE/dx energy loss of a fast particle traversing a fully ionized, 
but non-relativistic dilute plasma and the Lamb shift of hydrogen-like atoms. We shall 
first work out the energy loss problem in some detail because it involves rather elementary 
physics, and a self-contained exposition can be presented within a short space. The plasma 
example also does share some common features with relativistic plasmas such as those that 
appear in QCD and weak interaction physics. Indeed, this simple problem serves as a 
useful test of the validity of methods used in those more complex problems. In this regard, 
it should be noted that although our non-relativistic, Abelian example is conceptionally 
simpler than those in the relativistic, non-Abelian gauge theory, the non-relativistic plasma 
involves a Debye length and a plasma frequency that are quite distinct parameters while, in 
the extreme relativistic theory, these parameters become essentially the same (in lowest-order 
perturbation theory.) Thus a verification of the treatment of different physical processes may 
be more difficult in the relativistic theory, and the non-relativistic case may well serve as a 
quite useful testing bed for methods used in the relativistic case. Here our aim is to explain 
our method, and although we shall re-derive results that have been obtained before, the 
derivations will clearly describe and illustrate the power of the method. It has been applied 
recently to plasma energy-loss problems in which the constants under the logarithms were 
not known 0, [@. 

Dimensional regularization is widely employed in relativistic quantum field theory to 
make otherwise divergent expressions finite and well behaved, and then to implement the 
renormalization procedure. It should be emphasized that here we are making use of a 
different and novel application of continuation to spatial dimensions z/ 7^ 3 to compute results 
that are always well-defined and finite at the physical z/ = 3 dimension.^] We are not using 
dimensional continuation to render infinities finite so as to perform renormalizations as one 
does in quantum field theory. Moreover, our purpose here is to introduce and describe this 
new application of dimensional continuation by two very different but well-known physical 
problems in order to illustrate its range of application. What is new is the method. 

It is worth first illustrating our method with a trivial mathematical example, the behavior 
of the modified Hankel function Ky{z) in the small argument z limit with the index u also 
small. The argument z will play the role of the small parameter in our work below; the 
index v will play the role of the dimensionality except that in this simple Bessel function 
example we shall examine the region where v is near zero, not three. In general, the Hankel 



In the Lamb shift example at the end of this paper, we work with renormalized quantities that 
are well-defined and finite. 
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function has the integral representation 



Although k is simply a dummy integration variable, it is convenient to think of it as a wave 
number or momentum variable. When z is small, exp | — | + ^)} may be replaced by 1 
except when one or the other of the factors exp{—z k/2} or exp{— 2;/(2 k)} is needed to make 
the k integration converge in the neighborhood of one of its end points. When u is slightly 
less than zero, the integral (|1|) is dominated by the small k, "infrared or long-distance", 
region. In this case, only the exp{— 2;/(2 A;)} factor is needed to provide convergence, and 
we have 

z/ < : 

The variable change k = z/{2t) places this integral in the form of the standard representation 
of the gamma function, and we thus find that the leading term for small z in the region 
z/ < is given by 

z/ < : 




where 7 = 0.5772 ■ ■ ■ is Euler's constant. Note that the second line describes the behavior 
for u < near u = including the correct finite constant as well as the singular pole term. 

When u is slightly greater than zero, the integral (|l]) is dominated by the large k, "ul- 
traviolet or short- distance" regions. In this case, only the exp{—z k/2} factor is needed to 
provide convergence, and we have 

z/ > : 

The integral again defines a gamma function, and so 
z/ > : 

^.W^^ g)"^(l-^7), (5) 

with again the result containing the correct finite constant as well as the singular pole term. 

The result (0) for z/ < can be analytically continued into the region z/ > 0. In this 
region it involves a higher power of z than that which appears in the other evaluation (^, 
and hence this analytic continuation of the leading result for z/ < into the region z/ > 
becomes sub-leading here. Similarly, the result for z/ > may be analytically continued 
into the region u < where it now becomes sub-leading. An examination of the defining 
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integral representation (P shows that these sub-leading analytic continuation terms are, in 
fact, the dominant, first-subleading terms.0 For z/ > one term is leading and the other 
subleading, while for u < their roles are interchanged. Thus their sum 

contains both the leading and the first subleading terms for both u > and u < 0. In 
the limit u ^ the ("infrared" and "ultraviolet") pole terms in this sum cancel, with the 
variation of the residues of the poles producing a logarithm, yielding the familiar small z 
result 

Ko{z) = -\n{z/2)-^. (7) 

It must be emphasized that the correct constant terms [In 2 — 7] are obtained by this method 
in addition to the logarithm — In 2; which is large for small z. The result was derived 
from the analytic continuation of results that were easy to compute in one or the another 
region where either "infrared" or "ultraviolet" terms dominated. This is the essence of our 
method. Of course, the general result (H) could be obtained by a more careful computation 
of both the leading and first-subleading terms in either of the separate i/ > or z/ < regions 
as was done in the previous footnote. In some of the physical examples that we shall work 
out, such an extraction of the subdominant terms would be very difficult indeed, although 
possible in principle. Thus our method acquires real power for the physical problems. 

In the physical examples described below, one could object that we do not explicitly prove 
that larger subleading terms are not present. In the limit of the typical small parameter 
z ^ that concerns us (with now z/ — 3), for our physical examples we shall find (after 
extracting some overall factor of z) leading terms of order z^^~'^^ for u > 3 and z^'^~^^ for 
u < 3. One may then ask if an additional term that has a power dependence between 
z^'^~^^ and z^^^"^ can appear. However, simple dimensional analysis shows that such terms 
of intermediate order cannot appear. The point is that these physical examples involve only 
two different physical mechanisms that dominate at large and small scales, and that the 
two different mechanisms involve different combinations of the basic physical parameters 
and hence give quite different dependencies on the small parameter when the dimension u 



^For example, subtracting the leading term (^) for < from the integral representation (|T|) gives 
K^iz) - i (I) ' r(-z.) = '^f^k^ [e-'"' - 1] e-^l^^^^ . 

For > > — 1, the integral on the right-hand-side of the equation converges when the the final 
exponential factor in the integrand is replaced by unity, the 2; — > limit of this factor. Hence this 
final factor may be omitted in the evaluation of the first sub-leading term. A partial integration 
presents the result as 

whose evaluation gives precisely the analytic continuation of the leading term (Isl) for 1/ > 0. 
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departs from z/ = 3. Incidentally, it should go without saying that the physical examples 
that we shall study have a basic theoretical description that is valid for a range of spatial 
dimensions v about z/ = 3. 

II. ENERGY LOSS IN A PLASMA 

The usual method for obtaining the energy loss for a charged particle moving through 
matter is to divide the calculation into two parts: The long-distance, soft collisions and the 
short- distance, hard collisions. Collective effects are important in the long-distance part, 
and it is evaluated from the j ■ E power loss of a particle moving in a dielectric medium. The 
hard collisions are described by Coulomb scattering. The rub is to join the disparate pieces 
together. For the case of classical scattering, this is often done by computing the energy loss 
in Coulomb scattering out to some impact parameter, and then adding the j ■ E energy loss 
for all larger impact parameters. Although such methods do yield the correct large logarithm 
without much difficulty, the logarithm of the ratio of the two scales which is large, the purely 
numerical constants (which one expects to be of order one) that accompany the logarithm 
are harder to compute. Here we describe an easily applied method that yields a unique 
result - the result including the constants in addition to the large logarithm. The new idea 
is to compute the energy loss from Coulomb scattering over all angles, but for dimensions 
> 3 where there are no infrared divergences. A separate calculation of the energy loss 
using the j ■ E heating is done for z/ < 3, where the volume integration may be extended 
down to the particle's position without encountering an ultraviolet divergence. Both of these 
results have a simple pole at z/ = 3, but they both may be analytically continued beyond 
their initial range of validity. In their original domain of dimension z/, both calculations are 
performed to the leading order in the plasma density. As will be seen, although the Coulomb 
scattering result is the leading order contribution for z/ > 3, it is of subleading order when 
z/ < 3. Conversely, the j • E heating is subleading for z/ > 3 but leading for z/ < 3. Hence, the 
sum of the two (analytically continued) processes gives the leading and (first) subleading 
terms in the plasma density for all dimensions z/, and thus, in the limit of this sum at z/ = 3, 
the pole terms must cancel with the remainder yielding the correct physical limit to leading 
order in the plasma density. 

The fully ionized, classical plasma with which we are concerned is described exactly 
by a coupled set of kinetic equations, the well-known BBGKY hierarchy as described, for 
example, in Section 3.5 of ref. 0. We are interested, however, in the computation to leading 
order in the plasma density of the energy loss of a fast particle traversing the plasma. The 
correct equations that govern the leading order low-density behavior change as the spatial 
dimensionality v changes. For z/ < 3, the long-distance, collective effects dominate, and 
the formula derived by Lenard and Balescu applies [^, |Q. This formula describes the 
interaction of the various species that the plasma may contain. In the limit in which one 
species is very dilute, as is our case in which we examine the motion of a single, fast "test 
particle" moving through the plasma, the energy lost in the particle motion is described 
by its j • E Joule heating with the background plasma response given by the permittivity 
of a collisionless plasma. On the other hand, when the spatial dimension v is greater than 
3, the short- distance, hard Coulomb collisions dominate. For these dimensions, the leading 
low density limit of the BBGKY hierarchy is described by the familiar Boltzmann equation. 
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The Boltzmann equation is derived, for example, also in Section 3.5 of ref. 0. We use the 
Boltzmann equation to obtain the leading order energy loss rate when u > 3. Again, since 
we are concerned with the motion of a single, fast "test particle" , the Boltzmann equation 
reduces to the product of the energy loss weighted cross section times the plasma density. 
The derivations that we have just described, which start from first principles, justify the 
methods outlined in the previous paragraph, the methods that we shall use. 

Since we are only interested in describing the new method, we simplify the discussion by 
treating only the electrons in a classical plasma (electron recoil gives the dominant energy 
loss since they are light), and by taking the moving projectile velocity Vp to be much larger 
than the electron velocities in the plasma so that the latter may be neglected relative to 
Vp. We shall assume, however, that the projectile velocity is small in comparison with the 
velocity of light so that this particle produces a simple Coulomb field (as modified by the 
plasma) and that non-relativistic mechanics applies. 



A. 1/ < 3 

We first compute the j ■ E heating with z/ < 3. The electric field E is produced by the 
point projectile of charge Cp moving with velocity Vp which gives the charge density 

p{r,t) =ep6^''\r-Vpt) 



In our non-relativistic limit, the electric field is curl free while V ■ D = p, and so 

p F.(k.'k.v,) {'(k ■ ■■ - k . V, i)} , (9) 

where e(k, uj) is the wavenumber and frequency dependent electric permittivity of the plasma. 
Note that we use rationalized Gaussian units so that electrostatic potential of a point charge 
in three dimensions has the form = e/(47rr). Since the current of the projectile is given 
by j(r, t) = Vpp(r, t), this energy loss mechanism gives dE/dt = —CpWp ■ E(vpt, t), or 

dE< _ 2 r jd^k) i k ■ Vp 



dt PJ {2ny k^e{k,k-Vp) 

The electric permittivity is the boundary value of an analytic function, e(k, u) = e(k, cj + ir^), 
1] 0+, with e(k, 2;) = e(k, — 2) an even function of z which also depends only on the 
modulus |k| of k (by rotational invariance). [This is demonstrated in the explicit form (|T2p 
below.] Thus, in view of the factor k ■ Vp which extracts the odd part of the rest of the 
integrand in Eq. (p!oD , we may write this energy loss as 

d_E,^ rjd^^ J -k-vp 1 
dt PJ {27ryk^ \e{k,k-Vp)j ' ^ ^ 

in which Im denotes the imaginary part. 
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For our calculation to leading order in the plasma density, the permittivity function may 
be taken in the first (one-loop) approximation (the ring graph of quantum statistical 
mechanics) 

e(ku) = l-^f ^ n (V) 2[(p + ^k)V(2m.)-pV(2m.)] 

^ ' ^ k^J {27rhy h\uj + ir]y-[{p + hky/{2me)~py{2me)]^' 

where the 77 — >• 0"*" in the denominator corresponds to a retarded response, and where 

(d^p) 



ne(p) = He (13) 



is the electron number density, and rrie is the mass of the electron. Since we are examining 
the limit in which the projectile is moving very rapidly with respect to the average thermal 
velocity of the electrons, the electron motion in the plasma may neglected. This corresponds 
to setting p = in the remainder of the integrand in Eq. (|12D , which gives 



<^^^) = ^-T-r-^^^Tn:^77^2^ (14) 



{io + ir]Y — {hk'^/2m^ 
where Ue is the plasma frequency for the electrons defined by 
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Tie 



e 

nip 



(15) 



With u = k ■ Vp, two length scales appear, Vp/ue and Ti/imeVp). It is easy to check that 
the scale for the wave number integration in the energy loss (|lTD (with, of course, v < 
3) is set by the former, classical length. Hence the latter quantum length appears as a 
correction involving the dimensionless parameter h:^uj1/ {meVpY' . Since ujI is proportional to 
the electron density, and we are working to leading order in this density, we must omit this 
small parameter and use the purely classical limitQ 



= (16) 

This is the limit to be used in the energy loss ([TT|). In this limit. 

Hence, performing the integration over the component of k parallel to Vp, and writing 
dx = Vpdt gives 



^Note that our dimensional regularization method yields the correct leading-order result, unac- 
companied by any higher-order terms that only give a part of the higher-order corrections and thus 
represent spurious corrections. 
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dx 2 J (27r)^-i cj2 + v^k^ ■ ^ ^ 

Exponentiating the denominator via 

D^^ = / dse-'^ , (19) 



interchanging integrals, performing the resulting z/ — 1 Gaussian k integrals, and recognizing 
the final s integral as a standard representation of the F function gives 



dE^ _el ( ul \ ""2' „ /3 - z/ 



dx 2 \AT[Vp ^ 



r , (20) 



or, with the neglect of terms which vanish when i/ — 3, 

(21) 



dE< _ elujl ( uol \ { I 7 



dx AnVp \AnVp J [S — v 2 



The pole in this expression, which becomes negative when u > 3, corresponds to the ultra- 
violet divergence which appears when i/ — >■ 3 in the wavenumber integral (piSl) . 



B. u>3 

We turn now to the u > 3 case where the energy loss is computed by single-particle 
scattering. By the conservation of energy, the energy loss in the scattering of the projectile 
velocity Vp on electrons whose initial velocity may be neglected is 



= — ^ Vj, -Vp = — f e , (22) 



me 

2 



where fg is the speed of the scattered electron. Since the initial electron has negligible 
momentum, this can be written in the invariant form AE = g^/(2me), where q is the 
electron momentum transfer in the scattering process. With the initial electron at rest, the 
differential rate of scattering is VpTieda, where Ue is the electron density in the plasma and 
da is the cross section element. Since dx = Vpdt, the energy loss for z/ > 3 is given by 



dEy Up 



dx 2mp 



daq\ (23) 



We first evaluate this scattering contribution when the interaction is weak, when rj = 
Cpc/hvp -C 1. In this case, the quantum- mechanical Born approximation result is appropriate 



with, in z/ > 3 dimensions. 



d^Bq'= I ^2r.n5{^-^\y^\ - q\ (24) 




{2'KhY \2m 

Here (1/m) = (l/mg) + (l/m^) defines the reduced mass m and v is the relative velocity 
between the electron and the projectile. Writing q^ = Am'^v^ sin^ 6'/2 , and 
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{d^p') = mp'^"-'^^ d{p'^/2m) fi^„2 9 d9 , (25) 

with sin'^^^ 6 = [2cos^^/2 sin6/2Y~^ , and noting that the sohd angle is given by 



27r rf^^ 



(26) 



we get, on setting x = 9/2, 



rfcr^ = F% -yfzTY / dx cos^"^ x sin^^ x • (27) 



27rt>2 V nh 



The integral which appears here has the value (z/ — 3)^^ + 0(z/ — 3) as one can show by 
dividing it into two parts with a suitable partial integration or by expressing it in terms of 
the standard integral representation of the Beta function. Since we neglect the motion of 
the initial electron, the relative velocity v may be replaced by the projectile velocity Vp, and 
so usmg the result in Eq. (|2|) gives 



dx Attv^ V vrr / U - 3 2 



p 



(28) 



The pole in this expression, which becomes negative when z/ < 3, corresponds to the infrared 
divergence of the momentum integral (|2^) in the z/ — > 3 limit. 



C. z/ = 3 

When the result (pS]) is added to that in Eq. ( pi]) the divergent pole terms cancel, and 
the physical limit u —>■ 3 is 

dx A-nvl \ hixJe ) ' 

For small r/, this is the correct result to leading order in the plasma density. Instead of using 
the plasma density for the demonstration, it is equivalent to use the linearly related plasma 
frequency tUg. We have computed the leading and subleading terms in this quantity. The 
result (^) for dE^/dx involves ujI x uj^'^~^^ while the result (28) for dE^/dx involves just 
tUg. Hence, for z/ < 3, (|21|) is leading and (|28| ) is subleading, while for z/ > 3, their roles are 
reversed. Thus, in either region the sum of the two contributions contains both the leading 
and (first) subleading terms, and so the limit of the sum at the physical dimension u = 3 
yields the correct result to leading order in the plasma density. 

The result (PP]), including the proper constants inside the logarithm, may also be es- 
sentially obtained by applying the j ■ E heating formula (|1^) directly in three dimensions 
with the use of the single-ring graph quantum form (|T^ of the dielectric function in the 
limit in which the electrons in the plasma are taken to have negligible velocity. Placing this 
function in Eq. (0) with z/ = 3, writing k ■ Vp = cos^ in the resulting delta function, and 
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using this delta function to eliminate the polar angle 6 of the solid angle integration gives 
a remaining integral over the magnitude k of J dk/k. The leading terms for small cjg of 



the upper and lower limits of this logarithmic integral give the result (pQ]), except that the 
correct reduced mass m in Eq. (^) is replaced by the electron mass rrie since the current j 
describes the motion of a very heavy projectile particle. This sort of calculation was done 
some time ago by Lindhard [0. Although the reduced mass correction is negligible when 
the projectile is a heavy ion, it does represent a conceptual shortcoming of the quantum- 
corrected, joule heating treatment. Moreover, this treatment completely breaks down when 
the projectile is itself an electron. This sort of dielectric treatment is also restricted to the 
case of a cold plasma whose electron velocities are much less than that of the projectile. On 
the other hand, our method is easily extended [|l|] to treat the case of a hot plasma where 
this restriction is not imposed, and again a complete calculation can be performed which 
includes the constants in addition to the logarithm. 

Although, as we have just seen, the j ■ E calculation can be improved to obtain the 
correct energy loss (except for the replacement of the reduced by the electron mass), with 
the computation always done in three dimensions, we do not know of a similar improvement 
of the Boltzmann equation in three dimensions which yields the correct result. One might be 
tempted to replace the Coulomb potential by the screened Debye potential. This alteration 
changes the factor in the cross section formula by 

(30) 

in which = e^n/T is the squared Debye wave number for the plasma. This alter- 
ation removes the long-distance infra-red divergence, and the cross section formula ( p^ now 
converges in three dimensions. Exponentiating this denominator using the integral repre- 
sentation ([T9|) with an additional factor of s in the integrand to produce the square makes 
the remainder of the calculation easy, and one finds that 




dEo al^l 



dx An Vp 



In 



2mVp \ 1 



hK J 2 



(31) 



Although the constant out in front of the logarithm is again the correct over-all constant, 
the argument of the logarithm is quite different from the correct form given in Eq. (p9|). 
This should have been expected at the outset because Debye screening describes the static 
screening of a particle at rest in the plasma, not a dynamical screening of a fast moving 
particle which is the case that we are examining. As far as I know, such a dynamical 
screening within a Boltzmann equation context cannot be done. 

Our method can be used to extend the result ( PD| ) to arbitrary values of r] = eCp/ {An hvp), 



always retaining the correct additional constants. To do this, we use some clever mathemat- 
ics of Lindhard and Sorensen but in a manner which justifies that these constants have 
been kept. Namely, we compute 



. dE^ n„ 
A- 



dx 2mp 



J {da - das) ■ (32) 



This difference is well behaved in the limit u ^ 3 since the pole at = 3 produced by the 
cross section integral comes from soft, infrared physics which is completely contained in the 
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Born approximation das- Hence the three-dimensional partial wave decomposition of the 
scattering amplitude may be used, and then standard manipulations yield 



{da - dae) = 2t:Ti^ + 1) 

1=0 



2 — g2i[(5;-<5(,+i)] _ g-2i[<5;-(5(,+i)] 
2 — g2«[<5i-<5(i+i)] _ g-2i[5i-<5(,+i)] 



B 



(33) 



For the Coulomb potential 



T{l + l+i7]) 

r(/ + l-zr/)* 



(34) 



where the phase is independent of /. Using T{z + 1) = zT{z) and a little algebra, we find 
that 



/oo 
[da - d(jB) = ^-nrfTi^ ^ 



+ 



1 



I + 1 + 17] l + l-ir] l + l 



f4_ 

Att 



2[Re^(l + z77) +7] 



(35) 



where iplz) is the logarithmic derivative of the gamma function, ip{z) = r'{z)/r{z), and Re 
denotes the real part. Recalling the definition (plSf ) of the plasma frequency, we now have 



dx 



47r 



[ReV^(l + ir/) +7] 



(36) 



with the energy loss for all r] values given by 



dE _ dEgm _^ ^ ^^> 



An Vp 



- [ReV'(l + ^r/) +7] 



In the classical case, r] = eCp/ {An hvp) becomes large. Using the limit 
1^1 ^ 00 : tpil + z) = \nz + 0{z-'^) , 
Eq. ( ^71) yields the classical form 

dEci 



dx 



A-KV^ \ 



(37) 



(38) 



(39) 



This result, including the proper constant 2e that appears within the logarithm, was 



obtained long ago by Kramers [|T0| . It may also be obtained directly ||TT| with our dimensional 
continuation methods by using the classical Coulomb scattering cross section for dimension 
z/ > 3 in the scattering energy loss expression ([23|). 
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III. LAMB SHIFT 



Essentially the same method applied here has been used before in my Quantum Field 
Theory book to calculate the Lamb shift for hydrogen-like atoms, with the small pa- 
rameter role of the plasma density replaced by the nuclear charge Ze. That exposition, 
however, was presented in a somewhat mystical manner, and it unfortunately did not bring 
out the essence of the method. This will be rectified now and the process will provide an- 
other example of how the method works. Section 8.7 of ref. explains in detail how a 
radiative energy correction may be expressed as a matrix element of the electron self-energy 
operatoiQ "^[E) in Coulomb, bound-state Dirac wave functions of energy E. The Lamb 
shift is an energy difference that has both infrared and ultraviolet contributions just as in 
the more elementary plasma energy loss example explained above. In v spatial dimensions, 
these become two distinct physical processes that scale in different ways with a characteristic 
atomic energyQ Removing a common overall factor, the ultraviolet contribution behaves 
as while the infrared contribution goes as . Since the energy E vanishes when Ze^ 
vanishes, we may take E (implicitly divided by some fixed energy scale to yield a dimen- 
sionless number) as our small parameter. Thus, just as in the previous plasma case, the 
infrared part dominates when i/ < 3, the ultraviolet part dominates when i/ > 3, the sum 
of the two contributions analytically extended in the vicinity of = 3 always contains both 
the dominant and leading sub-dominant terms, and so the — 3 limit of this sum yields 
the correct leading-order Lamb shift. To tame the infrared divergences which are prevalent 
when z/ < 3, the binding of the electron must be accounted for. To tame the ultraviolet 
divergences that may appear when i/ > 3, a relativistic treatment must be made. We turn 
now to sketch this calculation. The detailed expressions that we shall need are derived and 



presented in ref. [12 



With soft (virtual) photon exchange, the leading terms are given by the non-relativistic 
limit of the electron motion. In this non-relativistic limit, the calculation is most easily 
performed in the radiation gauge. The Coulomb self-energy contribution is removed by 
a mass renormalization. The photon exchange contribution to the electron self-energy is 
properly renormalized by a subtraction so that it vanishes for a free particle. An elementary 
computation yields 



^In this formulation, the vacuum polarization contribution to the Lamb shift appears as a separate 
modification of the Coulomb potential. This is a simple correction which does not involve an 
interplay between long and short distances that concerns us in this paper, so we omit the effect of 
vacuum polarization here. 

^Simple dimensional analysis shows that characteristic atomic energy in v spatial dimensions is 
given by f = {fi^ /m) {Z'^e^ m? /fi^)^^~^^ ^ , which reduces to the familiar scale £ = Z'^e^ rajh} in 
three dimensions. Since the scaling behavior of S(£') is more simply expressed in terms of E rather 
than a dimensionless parameter formed from Ze^, we use E as our small parameter. 
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FIG. 1. Feynman diagram for the electron self-energy operator of Eq. ([40D. The waving line 
stands for the transverse virtual photon exchange of the radiation gauge. The narrowly spaced 
straight lines stand for the non-relativistic electron (operator) propagator in the nuclear Coulomb 
field. 



E<(i?) 



H-E 



P 



2k^ mc H — E + khc — ie mc 



(40) 



Here H is the non-relativistic Hamiltonian for the hydrogen-like atom with nuclear charge 
Ze. This result, which may be obtained from old-fashioned second-order time- dependent 
perturbation theory, involves the atomic Coulomb exchange to all orders as shown in Fig. 1. 
It is just Eq. (5) of the original Lamb shift paper of Bethe [|l^ except that it is written 
in u rather than 3 dimensions. It is also essentially Eq. (8.7.43) of ref. ||12|. Performing 



the integrations (as are explicitly done in the ref. |]T2|) gives, with the neglect of terms that 
vanish at u = 3, 
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mc 



(41) 



In describing the scaling of the results with respect to the small parameter S, we implicitly 
consider matrix elements of the self-energy operator in a bound-state energy eigenfunction 
and omit the scale associated with the two p operators that always appear in the expressions. 
Here, the two p flank the operator {H — Ey~^ which has the characteristic atomic size S"'"^. 
The operator {H — Ey~^ has this typical scale for any intermediate state when a complete 
set of intermediate states are inserted within the matrix element. Since the whole expression 
converges, it has the size S'^~'^. Thus we confirm the that the leading term for z/ < 3 in the 
Lamb shift is of order S'^~'^ as stated before. The divergence that appears when u approaches 
3 is, in view of the structure of the integral (^), an ultraviolet divergence. 



B. u>3 

In leading order, for spatial dimensionality greater than three, the leading correction 
to the electron self energy involves a single Coulomb interaction between the electron and 
the nucleus as shown in Fig. 2. The correction entails the radiative modification of this 
interaction which are described by the order term in the electric form factor -Fi(g^) of 
the electron, E[{0) (the charge radius), and a relativistic anomalous magnetic moment effect 
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A 

FIG. 2. Feynman diagram for the electron self-energy operator corrected by a single Coulomb 
exchange with the nucleus. The waving line stands for the transverse virtual photon exchange of 
the radiation gauge. The dashed line ending in a cross stands for the single Coulomb exchange 
with the nucleus. 



given by the magnetic form factor at zero momentum transfer, F2{0). The wave number 
factor gives the Laplacian of the Coulomb potential, V. Using the Dirac equation to 
reduce the -^2(0) contribution yields a term involving V'^V and a spin-orbit interaction. The 
Schrodinger energy eigenstate matrix elements of the operator h V are the same as those 
of the operator 2 p ■ {H — E) p. The results are derived in ref. |]12[ and presented there in 
Eq's. (8.7.31) and (8.7.36), which we now write together as 



S>(i?) 



Stt Att he 



-3 2 



mc 
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mc mc 



Ait he Airm'^c^ 



mc 



u-3 



a ■ i[p,V] X p, 



(42) 



where again terms that vanish at u = 3 are omitted. Note that, with our conventions, this 
result scales as S, again confirming an assertion made before. The divergence in the first 
term on the right-hand side of Eq. (^2]) when u approaches 3 comes from the contribution 



of F{(0) which contains an infrared divergence in three spatial dimensions |T2 



C. i/ = 3 

Since the two effects that we have listed contain the leading and first sub-leading terms 
for u near 3, their sum 

E(E) = E<(E) + S>(E) (43) 

evaluated at u = 3 must produce the Lamb shift to leading order. Indeed, the pole terms 
cancel as they must with the variation of the residues producing a logarithmic contribution, 
and one finds that 
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S(E) = |^ ^.{H-E)L 
Sir mc [ 



mc 



2{H -E 




o ■ hVV X p , (44) 



where now a = e^/ {Air he) ~ 1/137 may now be identified with the fine structure constant. 
This is the famihar formQ of the Lamb shift operator. This form appears in Eq. (8.7.63) 
in Ref. |jl2|, and its consequences are explained there. The correct factor of 5/6 has an 
interesting history in the computation of the Lamb shift, as related in footnote 13 of ref. 



It is worth noting that this calculation of the Lamb shift using our dimensional continu- 
ation method is simpler than that done using the methods of effective quantum field theory 



151 which would entail an additional matching calculation. It is much simpler than the con- 



ventional, old-fashioned method which utilizes a fictitious photon mass and a cumbersome 
joining process with an intermediate, non-covariant photon momentum cutoff.Q 
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^The imaginary part gives the width or lifetime of the level. 



^This hoary procedure is still presented in detail in modern texts on quantum field theory. See, 
for example, Section 7-3-2 of Itzykson and Zuber |16], or Section 14.3 of Weinberg [p^]. 
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